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We employ a bispinor gap equation to study superfluidity 
at nonzero chemical potential: jj, ^ 0, in two- and three-colour 
QCD. The two-colour theory, QC2D, is an excellent exemplar: 
the order of truncation of the quark-quark scattering kernel: 
K, has no qualitative impact, which allows a straightforward 
elucidation of the effects of fi when the coupling is strong. In 
rainbow-ladder truncation, diquark bound states appear in 
the spectrum of the three-colour theory, a defect that is elim- 
inated by an improvement of K. The corrected gap equation 
describes a superfluid phase that is semi-quantitatively simi- 
lar to that obtained using the rainbow truncation. A model 
study suggests that the width of the superfluid gap and the 
transition point in QC2D provide reliable quantitative esti- 
mates of those quantities in QCD. 



I. INTRODUCTION 

In the application of Dyson- Schwinger equations |l]] 
(DSEs) extensive use has been made of models based on 
the rainbow-ladder truncation, with contemporary vari- 
ants Q providing improved links with QCD. This trun- 
cation is also implicit in the class of model field theories 
with four-fermion interactions, such as the Nambu-Jona- 
Lasinio (NJL) model f| and the Global Colour Model @, 
which have been used successfully in describing aspects 
of the strong interaction. Such models admit the con- 
struction |ij of a meson-diquark auxiliary-field effective- 
action, which is important in developing an understand- 
ing of nucleons using the relativistic Fadde'ev equation 
|@|. In addition, it is immediately apparent that the ac- 
tion's steepest-descent equations admit the possibility of 
diquark condensation; i.e., quark-quark Cooper pairing, 
and that was first explored using a simple version of the 
NJL model @. 

A nonzero chemical potential: \x ^ 0, promotes Cooper 
pairing in fermion systems, and earlier and independent 
of these developments in QCD phenomenology, the pos- 
sibility that it is exhibited in quark matter was consid- 
ered m using the rainbow-ladder truncation of the gap 
equation. A quark-quark Cooper pair is a composite 
boson with both electric and colour charge, and hence 
superfluidity in quark matter entails superconductivity 
and colour superconductivity. However, the last feature 
makes it difficult to identify an order parameter that 
can characterise a transition to the superfluid phase: the 
Cooper pair is gauge dependent and an order parameter 
is ideally describable by a gauge-invariant operator. 

Determining the (T, /x) phase diagram of QCD is an 
important goal. At (T, fx) — there is a quark-antiquark 



condensate: (qq) ^ 0, but it is undermined by increasing 
T and [i, and there is a domain of the (T, /x)-plane for 
which (qq) = 0. Increasing T also opposes Cooper pair- 
ing. However, since increasing fi promotes it, there may 
be a (low-T,large-/i)-subdomain in which quark matter 
exists in a superfluid phase. That domain may not be ac- 
cessible at RHIC, which will concentrate on fj, ~ where 
all studies indicate that QCD with two light flavours ex- 
hibits a chiral symmetry restoring transition: (qq) — > 0, 
at T ~ 150 MeV. However, it may be discernible in the 
core of dense astrophysical objects which could un- 
dergo a transition to superfluid quark matter as they 
cool, and in baryon-density-rich heavy ion collisions at 
the BNL-AGS and CERN-SpS §]. An exploration of this 
possibility using numerical simulations of lattice-QCD is 
inhibited by the absence of: (i) a gauge-independent or- 
der parameter for the superfluid phase; and (ii) a sat- 
isfactory procedure for the numerical estimation of an 
integral with a complex measure, such as the fJ,j^0 QCD 
partition function. Consequently all the information we 
have comes from models. 

The rainbow-ladder truncation has the feature and de- 
fect that it generates a quark-quark scattering kernel, K, 
that is purely attractive in the colour antitriplet chan- 
nel, 3 C . It therefore not only yields a 3 C scalar diquark 
condensate but also 3 C diquark bound states ju]]; i.e., 
hitherto unobserved coloured quark-quark bound states 
with masses (in GeV) [ pd| : 



ud 



=0+ 



0.74, 



ud. 



0.95. 



m u n d 



1.5 



(1) 



(us = ds diquarks are also bound; e.g., TOq+ = 0.88. 
Colour-sextet bound states do not exist because K is 
purely repulsive in this channel, even in rainbow-ladder 
truncation ]i~0[].) All models employed to date in the 
analysis of quark matter superfluidity have this defect 
|l2j , and we are primarily concerned with the question of 
whether any model or truncation with such a flaw can be 
a reliable tool for exploring superfluidity in quark matter. 
In addressing this issue, it is important to compare QC2D 
with QCD because the same mechanism that provides for 
the absence of diquark bound states in the latter must 
guarantee their existence in QC2D, where the diquark is 
the baryon of the theory. In fact, it must ensure that 
fiavour-nonsinglet J p= ^ mesons are degenerate with J 
diquarks |l~3| ]. 

In Sect. II we describe a bispinor DSE (gap equation) 
that is particularly useful for studying quark and diquark 
condensation and, in Sect. Ill, employ it in the general 
analysis of QC2D and also to obtain quantitative results 
from a pedagogical model. In Sect. IV we focus on QCD, 
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and employ the model's analogue to exemplify the gap 
equation and its solution in rainbow truncation, and also 
when a 1 /^-suppressed dressed- ladder vertex correction 
is included. We summarise and conclude in Sect. V. 



II. A GAP EQUATION 

A direct means of determining whether a SU c (iV) 
gauge theory supports + diquark condensation is to 
study the gap equation satisfied byn 



V(p,n) := 
where, with u> 



(2) 



D(p,fi) 

-\*{ P ,-n) l5 \\ 



[ij] =Pi + in, 



A\p^) lb X\ 



(3) 



ij-pA(p 2 , JL ) + «74 u> w C (p 2 , JL ) + B(p 2 , Jf , ) 



{A 



1 



N C {N C 



l)/2} are the anti- 
symmetric generators of SU c (iV c ), and C = 7274 is the 
charge conjugation matrix: OyJC^ = — 7^ ; [C, 75] = 0. 
Using the gap equation to study superfluidity makes un- 
necessary a truncated bosonisation, which in all but the 
simplest models is a procedure difficult to improve sys- 
tematically. 

In addition to the usual colour, Dirac and isospin in- 
dices carried by the elements of T>(p, //), the explicit ma- 
trix structure in Eq. (|^) exhibits the quark bispinor index 
and is made with reference to 



Q(x) 



q(x) 
q(x) :=T 2 Cq T 

Q(x) :=(q(x) q(x) := q T C r 2 ) 



(4) 
(5) 



where {rl 



1,2,3} are Pauli matrices that act on 



the isospin index. Herein we only consider two-flavour 
theories, SU f(Nf — 2), because Nf does not affect the 
question at the core of our study, and focus on T = 0, 
since nonzero T can only act to eliminate a condensate. 
A nonzero quark condensate: (qq) 7^ 0, is represented in 



the solution of the gap equation by B(p' 



) ^ while 



diquark condensation is characterised by A l (p,/j,) ^ 0, 
for at least one i. 

The bispinor DSE can be written in the form 



A)(p,m) 



(6) 



Sn(p,M) S 12 (p,/x) 
7 4 Si 2 (-P,m)74 CSn(-p, ^) T C t 



1 In our Euclidean formulation: {7^, 7^} = 2<5 M „, 7^ = 7 M , 
p-q = Z)*=iP*9*> and tr D[757fi7i'7p7<T] = -4e M „ p<T , £1234 = 1- 



where in the absence of a diquark source term 
T> (p, n) = (i7 • p + m)rR - , 



(7) 



with m the current-quark mass. Here we have introduced 
additional Pauli matrices: {tq, a = 0, 1, 2, 3} with r° = 
diag(l, 1), that act on the bispinor indices. The structure 
of T,ij(p, n) specifies the theory and, in practice, also the 
approximation or truncation of it. 



III. TWO COLOURS 



As an important and instructive first example we con- 
sider QC2D. In this special case A l A^ = At 2 in Eq. (j|) 
and it is useful to employ a modified bispinor 



q{x) 

= t 2 q(x) 



(8) 



with Q2 the obvious analogue of Eq. (||), so that the La- 
grangian's fermion-gauge-boson interaction term is sim- 
ply Q2{x) ^gi^TQ Q%(x) A*(x) because SU C (2) is pseu- 

doreal; i.e., t 2 (— t c ) t t 2 = f c , and the fundamental and 
conjugate representations are equivalent. 

The gap equation at arbitrary order in the systematic, 
Ward-Takahashi identity preserving truncation scheme of 
Rcf. JTJ] is readily derived. For n = 0: C = A in Eq. (g), 
all the functions in the dressed-bispinor propagator are 
real and the rainbow truncation yields the gap equation 



2>(p) = 



17 • p + m + 



(9) 



^g 2 D, v {p~q)l/-f «%)7,y . 



(2t) 



Renormalisation is straightforward, even at [i ^ p5[ , 
however, since it is not relevant to our central theme, we 
neglect it here. To solve Eq. (^|) we consider a generali- 
sation (TJ of Eq. (|) 



V(p) = £7 • p A{p 2 ) + V(-tt) M(p 2 ) : 



(10) 



V(7r)=cxp| i75 £>V(p 2 )j =V(-tt)- 1 , (11) 

where {T 1 ' 2 ' 3 = t% <8 f}, T 4 = ® t], T 5 = r% <g> r^}, 
{T\Ti} = 2S ij , so that 



S( P ) = 



~ij-pA(p 2 ) + V{n)M{ P 2 ) 

p 2 A 2 (p 2 )+M 2 { P 2 ) 
-ry ■ pa v {p 2 ) + V(tt) a s (p 2 ) ■ 



(12) 
(13) 



(tt = (0, 0, 0, 0, — jn) produces an inverse bispinor prop- 
agator with the simple form in Eq. (0).) 

For ir(p 2 ) = constant, substituting Eq. ( JTo| ) into 
Eq. (|), yields 
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1 -p[A(p 2 )-l] 



T0f 9 2 D I1U (p - q) 7 M y 7 • q a v (q 2 ) 7^ y 



7W(p 2 ) -mV(Tr) = 



(14) 



(15) 



J ^g 2 D^(p~q) 1 ^a s (q 2 )^ 



It is clear from these equations that the gap equation 
in rainbow truncation is independent of tt in the chiral 
limit. As this result is true order-by-order in the trun- 
cation scheme of Ref . [[lij , it is a general property of the 
complete QC2D gap equation. Hence, if the interaction 
is strong enough to generate a mass gap, then that gap 
describes a five-parameter continuum of degenerate con- 
densates: 



(Q 2 v( 7r )g 2 )^o. 



(16) 



and there are 5 associated Goldstone bosons: 3 pions, 
a diquark and an anti-diquark, which is a well-known 
consequence of the Pauli-Giirsey symmetry of QC2D. 

For m ^ 0, it is clear from Eq. (15) that the gap 
equation requires tr/Q [T l V] = 0, i.e., in this case only 
(Q2Q2) 7^ 0. The Goldstone bosons are now massive but 
remain degenerate. 

The Landau gauge dressed-gauge-boson propagator is 



T{k 2 



(17) 



and to explore /1 ^ we employ a pedagogical model for 
the vacuum polarisation in QC2D: 



T 2 (k 2 ) = ^f, 2 5\k). 



(18) 



This form was introduced [0 for the modelling of con- 
finement in QCD. However, it is also appropriate here 
because the string tension in QC2D is nonzero, and that 
is represented implicitly in Eq. (18) via the mass-scale fj. 



Further, a simple extension of the model has been used 
efficaciously jl8|,[l9| as a heuristic tool for the analysis of 
QCD at nonzcn>[T, /1). The mass of the model's J = 1 
composites is a useful reference scale and for m = in 
rainbow-ladder truncation 



J=i 



(19) 



mj—i is weakly dependent on m, changing by < 2% on 
m G [0, 0.01] fj, while adding l/A^-suppressed vertex cor- 
rections produces an increase of < 10% Q. 

We now consider the [i ^ gap equation and suppose 
a solution of the form 



D (P,(J>) 75_A(p,/i) 
-j 5 A*(p,fi) D{p,n) 



(20) 



with D(p, /i) defined in Eq. (||) and D(jp,fi) :— 
C D(—p, n) CK In the absence of a diquark condensate; 
i.e., for A = 0, 



[U B (a),V(p,n)]=0, U B {a) := c 1 ' 



(21) 



which is a manifestation of baryon number conservation 
in QC 2 D. 

The inverse, S (p, fi) , is sufficiently complicated that 
it provides little insight directly. However, that can be 
obtained using Eq. (|l8|), which yields an algebraic gap 
equation. Using the rainbow truncation we find, at p 2 — 
\p\ 2 +p 2 4 = 0: 



A - 1 = \ fj 2 K {A (B* 2 - C* V) + A* |A| 2 } 
fi(C - 1) = f r) 2 K {C (B* 2 - C* V) - C* |A| 2 } 
B - m = t) 2 K {B (B* 2 - C*V) + B* |A| 2 } , 
A = fj 2 K {A (\B\ 2 + |C| 2 m 2 ) + A |A| 2 } , 



(22) 
(23) 
(24) 
(25) 

IAI 4 . 



within 1 = \B 2 -C 2 ^ 2 \ 2 + 2\A\ 2 (\B\ 2 + \C\ 2 ^ 2 ) 
These equations illustrate: the B <-> A degeneracy de- 
scribed above for (m, //) = 0; that A is real for all fi; and 
also the action of fi: enhancing the coupling in the A 
equation but suppressing it for B, which is how increas- 
ing (j, promotes diquark condensation at the expense of 
the quark condensate. 

For (m, u) = the rainbow gap equation is Eq. (JsJ) and 
with Eq. ( |l8[) the solution is: 



M 2 ( P 2 ) 
B 2 (p 2 )- 
A( P 2 ) = C( P 2 



& 2 (P 2 ) = \ I 



1 



Mi 



- Ap 2 , p 2 < \ 
otherwise 

P 2 <^ 
otherwise . 



(26) 



(27) 



The dynamically generated mass function, M(p 2 ), is tied 
to the existence of quark and/or diquark condensates and 
breaks chiral symmetry. Further, in combination with the 
momentum-dependent vector self energy, it ensures that 
the quark propagator does not have a Lehmann repre- 
sentation and hence can be interpreted as describing a 
confined quark |l],^(| . The interplay between the scalar 
and vector self energies is the key to realising confine- 
ment in this way pj[ , and is precluded in studies that 
discard the vector self energy. 

For fj, ^ and arbitrary p we solve the rainbow gap 
equation numerically, and determine whether quark or 
diquark condensation is stable by evaluating 



5P := P(fi, S[B = 0, A]) - P(ji, S[B, A = 0] 



(28) 



where the pressure is calculated using a steepest descent 
approximation p^] : 



P[S] = -TrLn[5] - iTr[(D - V )S] 



(29) 



SP > indicates that diquark condensation is favoured. 

The calculation of SP is facilitated by employing the 
^-dependent "bag constants" Jlffl 



3 



P(fi,S[B,A 



(30) 



0])-P(/x,S[P = 0,A = 0]), 



with B&(p) an obvious analogue. They measure the sta- 
bility of a quark- or diquark-condensed vacuum relative 
to that with chiral symmetry realised in the Wigner-Weyl 
mode. The (to, /i) = degeneracy of the quark and di- 
quark condensates is manifest in 



B B (0) = Sa(0) = (0.092 t)) 4 = (0.13 m J=1 ) 4 



(31) 



Increasing p, at m — and excluding diquark conden- 
sation one finds pq| chiral symmetry restoration at 



M 2c = 0.28?7 



(32) 

when Sb(/x) = 0; i.e., the pressure of the Wigner and 
quark-condensed phases are equal. However, for all /x > 

0: 5P > 0, B a (m) > °> with B A (Mfc A= °) = 
(0.20 mj = i) 4 > £>a(0). Therefore the vacuum is unsta- 
ble with respect to diquark condensation for all fi > 
and one has confinement and dynamical chiral symmetry 
breaking to arbitrarily large values. Of course, we have 
ignored the possibility that r) is /x-dependent. In a more 
realistic model, the /x-dependence of fj would be signifi- 
cant in the vicinity of M^' _ °> with fj — > as /x — > 00, 
which would ensure deconfinement and chiral symmetry 
restoration at large- fx. 

A ^ in Eq. @ corresponds to tt = (0, 0, 0, 0, |tt) 
in Eq. (|l6|); i.e., (Q^sTqC^) 7^ 0, and although the 
fi ^ theory is invariant under Q2 — > Ub(oc) Q 2 , Q 2 — > 
QiVb{— a), which is associated with baryon number 
conservation, the diquark condensate breaks this sym- 
metry: 

(Ga^Ga) (33) 

-> cos(2a) (Q 2 i75TQQ 2 ) - sin(2a) (Q 2 «75Tq<3 2 ) ■ 

Hence, for (to = 0, /x 7^ 0), one Goldstone mode remains. 

For to 7^ and small values of /x, the gap equation 
only admits a solution with A = 0; i.e., diquark conden- 
sation is blocked. However, with increasing fx a diquark 
condensate is generated; e.g., we find the following min- 
imum chemical potentials for diquark condensation 

to = 0.013 mj—i 
to = 0.13 mj—x 



M A ^o = 0.051 to j= i. 
fx^o = 0.092 TO, /= i. 



(34) 



Improving on rainbow-ladder truncation may yield 
quantitative changes of < 20% in the illustrative results 
provided by our model of QC 2 D. However, the pseudore- 
ality of QC 2 D and the equal dimension of the colour and 
bispinor spaces, which under ly the theory's Pauli-Giirsey 
symmetry, ensure that the entire discussion remains qual- 
itatively unchanged. QCD, however, has two significant 
differences: the dimension of the colour space is greater 
than that of the bispinor space and the fundamental 
and conjugate representations of the gauge group are not 
equivalent. The latter is of obvious importance because 
it entails that the quark-quark and quark-antiquark scat- 
tering matrices are qualitatively different. 



IV. THREE COLOURS 

In canvassing superfluidity in QCD we choose A l X\ 
A A 2 in Eq. (§) so that 



(35) 



D\\(jp,ii)I)\+D ± (j>,n)P ± 


A(p, /x) 75 A 2 


-A(p, -A*)7 5 A 2 


D n (p,fj,)P\\+D ± (p,^)P ± 



where P|| = (A 2 ) 2 , Pl+P|| = diag(l, 1, 1), and Dy ...D± 
are defined via obvious generalisations of Eqs. (||) and (||). 
The evident, demarcated block structure makes explicit 
the bispinor index. Here each block is a 3 x 3 colour 
matrix and the subscripts: ||, _L, indicate whether or not 
the subspace is accessible via A 2 . The bispinors associ- 
ated with this representation are given in Eqs. (^) and 
(||) , and in this case the Lagrangian's quark-gluon inter- 
action term is Q(x)igT p Q(x)A p \(x), 














(36) 



It is again straightforward to derive the gap equation at 
arbitrary order in the truncation scheme of Ref. p"4j ] and 
it is important to note that because 

D||(p,/x)P||+ J D x (p, / x)P x = A {|£>||(p,/x) (37) 
+lD ± (p, /x)} + ^A 8 {D|| (p, /x) - D ± (p, /x)} 

the interaction: T p [S(p, /x)r°, necessarily couples the ||- 
and _L-components. That interplay is discarded in mod- 
els that ignore the vector self energy of quarks, which 
is a necessary and qualitatively important feature of 
QCD ||jl9|gH|. 



A. Rainbow truncation 

Diquark condensation at /x = was studied in Ref. 
using a minor quantitative adjustment of the confining 
model gluon propagator defined via Eq. (|l§|): 

T{k 2 ) = 4ttV<5 4 (£0, (38) 

with which the rainbow-truncation gap equation is 

V(p, fj) = V (p, M ) + y^ry 2 Y a p S(p, (i) T a p . (39) 

Solving this and the ladder-truncation Bethe-Salpeter 



equation one obtains JL8 19 



hn 2 

2 '/ : 



(40) 
(41) 
(42) 



(qq) = (0.11 rjf. 
B B (p = 0) = (O.IO77) 4 , 



and momentum-dependent vector self energies, Eq. (|27 
which lead to an interaction between the ||- and _L- 
components of V that blocks diquark condensation. This 
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FIG. 1. /x-dependent "bag-constants" in the QCD model 
defined via Eq. Rainbow truncated gap equation - 

Dotted line: Bs(/x), short-dashed line: Ba(^)- At the in- 
tersection, where the system flips to the superfluid phase, 



4 ) = (0.75) Ss(O). Vertex-corrected gap equation 



- Solid line: Bb((J,), long-dashed line with circles: Ba{/j.)- At 
the intersection: Ba = (0.96) 4 Sb(0). The structure evident 
in 23a (m) 1s an artefact characteristic of Eq. 



is in spite of the fact that A°A 2 (— A a ) T = ±A a A a , which 
entails that the ladder-truncation quark-quark scattering 
kernel is purely attractive and strong enough to produce 
diquark bound states jl(J]. 

For \x ^ and in the absence of diquark condensa- 
tion, the model defined via Eq. ( |3g| ) exhibits jl8| coinci- 
dent, first order chiral symmetry restoring and deconfin- 
ing transitions at 



s,A=o 

t^c, rainbo 



= 0.28 77, 



(43) 



which is where Bb = 0. However, for /i ^ Eq. ( |39| ) 
admits a solution with A(p, /j) ^ and B(p, /j) = 0. 
SP in Eq. ( |2"8| ) again determines whether the quark- 
condensed or superfluid phase is the stable ground state. 
With increasing /i, Bb(^) decreases, very slowly at first, 
and Ba(i^) increases rapidly from zero. As illustrated in 
Fig. [I], that evolution continues until 

0.25 ?? = 0.89 /xf^bow- (44) 



B=0,A 
^c, rainbow 



where £>a (m) becomes greater-than Bb (/j). This signals 
a first order transition to the superfluid ground state^j 
and at the boundary 

(Q*75^A 2 Q) m= ^ : o,a^ = (0.65) 3 (QQ) M=0 . (45) 

These results are typical p5[ of rainbow truncation mod- 
els in which the parameters in the dressed-gluon propa- 
gator are tuned to yield the correct ir-p mass splitting. 




FIG. 2. Dashed line: A(z, /i^' A ) obtained in rainbow trun- 
cation with the QCD model defined via Eq. (^), plotted for 
a — as a function of p, where z = p (0, 0, sin a, ifi + cos a). 
As (i increases, the peak position shifts to larger values of p 
and the peak height increases. Solid line: A(z,/i = 0) ob- 
tained as the solution of Eq. (fl6|), the vertex-corrected gap 
equation, also with a — 0. 

The solution of the rainbow gap equation: A(p, fif' A ), 
which is real and characterises the diquark gap, is plotted 
in Fig. ^. It vanishes at p 2 — as a consequence of 

the | ] — I coupling that blocked diquark condensation at 

/i — 0, and also at large p 2 , which is the manifestation of 
asymptotic freedom in the model. 



B. Vertex corrected gap equation 

The next-order term in the gap equation corresponds 
to adding a 1 /A^-suppressed dressed-ladder correction to 
the quark-gluon vertex, and using Eq. ( |38| ) this yields 

(46) 



v(p, m) - x>o(p, m) + tev 2 s( P , M ) r; 
- 2§6^ 4 r p S(p, m) r» 5(p, M ) r; 5(p, M ) r* 

which is illustrated in Fig. [|. The kernel of the Bethe- 
Salpeter equation receives three additional contributions 
at this order. Their net effect is repulsive at timelike to- 
tal momentum and hence they prevent the formation of 
diquark bound states | p^,p6[ . The rf term in Eq. ( |4^ ) 
means that an algebraic solution cannot be obtained, 
however, a numerical solution is possible. For simplic- 
ity we only consider m = since m/0 opposes diquark 
condensation, as we saw in Sect. [II. At this order there 



is a A ^ solution even for fi = 0, which is illustrated in 




2 With 77 independent of /x, quark confinement, expressed as 
the absence of a Lehmann representation for S, persists in the 
superfluid domain. 



FIG. 3. Illustration of the dressed-ladder vertex-corrected 
gap equation, Eq. (|4q) . Each bispinor quark-gluon vertex is 
bare, given by Eq. pq). 
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Fig. ||, and 

m^(l.l) 2 m^ laddcr (47) 

(QQ) = (l.O) 3 (QQ) rainbow (48) 

B B = (l.l) 4 2?Jj ainbow , (49) 

where the rainbow-ladder results are given in Eqs. (f40|)- 
([42|), and 

(<3*75t|>A 2 Q) = (0.48) 3 (QQ) , (50) 

B A = (0.42) 4 £ B . (51) 

Unsurprisingly the quark-condensed phase is favoured. 
Precluding diquark condensation, the model exhibits co- 
incident, first order chiral symmetry restoring and decon- 
finement transitions^ at 

Mf' A =° = 0.77Mf r l= b ° ow - (52) 

Our numerical results^] for the /i-dependence of the 
"bag constants" are depicted in Fig. [I], which shows there 
is a transition to the superfluid phase at 

M f= > A =0.63 M f' A=0 , (53) 

and at the boundary (cf. Eq. ([45])) 

(Qi75^A 2 Q) M=0 . 63M f.A=o = (0.51) 3 (QQ> P=0 ■ (54) 

The ratio of the condensates increases by < 7% on 
/i 6 [0, fif =a ' A }. Quantitatively, the next-order correc- 
tion leads to a reduction in the critical chemical potential 
for the transition to superfluid quark matter but doesn't 
much affect the width of the gap. Qualitatively, the tran- 
sition occurs despite there being insufficient binding at 
this order to produce diquark bound states. 

Further insight is provided by solving the inhomoge- 
ncous Bcthc-Salpeter equation for the + diquark vertex 
in the quark-condensed phase. At fi = and zero to- 
tal momentum: P = 0, the additional contributions to 
the quark-quark scattering kernel generate an enhance- 
ment in the magnitude of the scalar functions in the 
Bethe-Salpeter amplitude. However, as P 2 evolves into 
the timelike region, the contributions become repulsive 
and block the formation of a diquark bound state. Con- 
versely, increasing /i at any given timelike-P 2 yields an 
enhancement in the magnitude of the scalar functions, 
and as (i — > /if ,A=0 that enhancement becomes large, 
which suggests the onset of an instability in the quark- 
condensed vacuum. 



V. EPILOGUE 

We have studied a confining model of QCD using a 
truncation of the Dyson- Schwinger equations that de- 
scribes well the ir-p mass splitting at (T,n) — and 
ensures that no diquark bound states appear in the spec- 
trum. Employing a criterion of maximal pressure, we ob- 
serve a first order transition to a chiral symmetry break- 
ing superfluid ground state, which occurs at a chemical 
potential approximately two-thirds that required to com- 
pletely eliminate the quark condensate in the absence of 
diquark condensation. Without fine-tuning, the super- 
fluid gap at the transition is large, approximately one-half 
of that characterising quark condensation. Thus, while 
completely changing the qualitative nature of the bound 
state spectrum in the model; i.e., eliminating unobserved 
coloured bound states, our vertex-corrected gap equa- 
tion yields a phase diagram that is semi-quantitatively 
the same as that obtained using the rainbow truncation. 
This bolsters our confidence in the foundation of current 
speculations |Q about the phases of high-density QCD. 

The procedure we used to improve the gap equation 
is equally applicable to two-colour QCD, which we anal- 
ysed with the help of a pedagogical model for the dressed- 
gauge boson propagator. Diquark bound states must ex- 
ist in QC2D because they are the baryons of the theory, 
and the truncation procedure ensures this, with the re- 
sult that flavour- nonsinglet J p= T mesons are degenerate 
with J diquarks. Using a straightforward, constructive 
approach, we saw that at fi = there are five Goldstone 
modes in QC2D, and that one of them survives at /J, 7^ 0. 
A nonzero current-quark mass opposes diquark conden- 
sation but for light-fcrmions there is always a value of the 
chemical potential at which a transition to the superfluid 
phase takes place. Our model studies indicate that in 
some respects; such as the transition point and magni- 
tude of the gap, the phase diagram of QC2D is quanti- 
tatively similar to that of QCD. This observation can be 
useful because the simplest superfluid order parameter is 
gauge invariant in QC2D and the fermion determinant is 
real and positive, which makes tractable the exploration 
of superfluidity in QC2D using numerical simulations of 
the lattice theory pq). The results of those studies can 
then be a reliable guide to features of QCD. 
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